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Abstract 

We investigate the total inclusive decay rate of the (ground state) Be meson within the 
framework of an operator product expansion in inverse powers of the heavy quark masses 
and subsequent matching onto nonrelativistic QCD. The expansion is organized as a series 
in the strong coupling and in powers of the heavy quark velocities in the B^, reflecting 
the nonrelativistic nature of a heavy-heavy bound state. In this aspect the character of 
the expansion differs from the more familiar case of heavy-light mesons. The framework 
incorporates systematically corrections to the leading b- and c-quark decays due to binding 
effects, as well as contributions from weak annihilation and Pauli interference. Based on this 
approach we find for the Be meson lifetime = (0.4 — 0.7) ps, the dominant mechanism 
being the decay of the charm constituent. 



PACS numbers: 13.25. Gv, 12.38.Lg, 11. 10. St 



I. INTRODUCTION 



The motivation for studying weak decays of heavy hadrons is essentially twofold. First, 
one aims at understanding basic properties of the weak interaction at a fundamental level, 
including the precise determination of CKM parameters. Second, systems containing heavy 
quarks allow us to test our understanding of QCD in an interesting limiting case where, due 
to the large mass scale involved, certain aspects of the dynamics simplify. Both topics are 
of course intimately connected as the analysis of weak decays of heavy hadrons is always 
faced with the problem of disentangling the interplay of the strong and weak forces. In this 
respect the Be meson, the lowest be quark bound state, is a particularly interesting system 
to study. Unlike in the case of heavy- light mesons, as for example B~^, Ed or i?^, the bound 
state dynamics of the Be can be systematically treated in a nonrelativistic expansion, which 
has proved very successful for the description of the cc and hb family. At the same time, and 
by contrast to the cc and bb ground states, B^. is stable against strong or electromagnetic 
decay due to its flavor content and disintegrates only via weak interactions. With these 
properties B^. is in fact a unique example, since the top quark lifetime is so short that the 
analogous £b or tc mesons do not exist. 

Several features of the be system have already been the subject of investigations in the 
past. A comprehensive analysis of the be spectroscopy and the strong and electromagnetic 
decays of the excited states has been given in Weak decay properties of the groundstate 
Bci semileptonic and various exclusive modes have also been discussed The lifetime of 

Bci tbci is briefly considered in where a rough estimate has been presented. An estimate 
of tb^ using a modified spectator model and information gained from the calculation of 
dominant exclusive modes is described in ||^. The wide range of lifetimes, r^^ = (0.4—1.2) ps, 
reported in these papers, reflects the various model assumptions on the modification of the 
free quark decay rates due to bound state effects. 

In the present article we discuss a systematic approach to computing the Be lifetime, 
which is based on the optical theorem for the inclusive decay rate, an operator product ex- 
pansion of the transition operator and a subsequent nonrelativistic expansion of the operator 
matrix elements in the B^. meson state. In the first step, the operator product expansion 
of the transition operator, we rely on the fact, that if both quarks can be considered as 
heavy (and their mass difference is also large), the energy release in the weak decay of either 
quark is large compared to the characteristic scale for the bound state dynamics. One may 
then expand in the ratio of these scales. Technically, this step copies the procedure for 
inclusive decays of heavy-light mesons, reviewed in |]^. In the second step, the expansion 
of matrix elements, we utilize that the bound state of two heavy quarks is nonrelativistic 
in first approximation. This approximation can be systematically improved by means of 
nonrelativistic QCD ^j, which organizes the evaluation of matrix elements in full QCD in 
an expansion in p/m^ and p/mc, where p = rribVb = rricVc ~ 1 GeV is the typical quark 
three momentum in the B^ meson. The finite set of matrix elements in nonrelativistic QCD, 
which incorporates all nonperturbative effects, must be determined from lattice calculations 
or, less rigorously, from potential models. At this point, our treatment differs from an anal- 
ogous one for heavy-light mesons, whose bound state dynamics is essentially different. In 
this case the scale relevant to the bound state is Aqcd, while in a heavy-heavy system an 
additional scale, p > Aqcd (for the case at hand), is dynamically generated. Consequently, 
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the importance of different operators is not ordered according to their dimension alone, but 
follows from the 'velocity scaling rules' derived in 0, when adapted to the case of two heavy 
quarks of different masses. 

The procedure outlined above when combined with the inevitable emission of hard gluons 
in the decay, results in a double expansion in as and vq = p/niQ, with both parameters small, 
if both quarks are sufficiently heavy. Although this feature is conceptually very attractive, 
it is not a priori clear whether in the realistic case the parameters will be small enough to 
guarantee a reasonable behavior of the expansion. First, in the c — > sud transition, the ratio 
of the typical quark momentum in the bound state and the energy release is not a very small 
number. Second, although the reduced mass for the he system falls in between the reduced 
masses of hh and cc, the c quark velocity in the Be meson is larger than in the J/^E', because 
the c quark has to balance the momentum of a heavier h quark. Thus relativistic corrections 
are expected to exceed those in the cc system. While it might appear that the meson is 
a rather marginal case for the operator product expansion, its convergence properties can 
only be properly assessed after an explicit numerical investigation. As we shall see, there 
are no obvious indications that the nonrelativistic expansion does not work. 

We estimate a short B^ lifetime (in comparison to earlier estimates), 

rs, = (0.4-0.7)ps, (1) 

the main uncertainty arising from the poorly known charm quark mass. Be mesons are 
expected to be produced at the LHC, if not before, in numbers that are sufficient to test this 
prediction. Recently the first B^. meson candidates have already been reported from LEP 
p!0| . The search for B^. at the Tevatron is summarized in [O . 



II. THE OPERATOR PRODUCT EXPANSION 

The optical theorem relates the total decay width of a particle to the imaginary part of 
its forward scattering amplitude. Applied to the Be meson total width Vb^ this relationship 
can be written as 

= —^(Be\T\Be) (2) 
where the transition operator T is defined by 

r = imi Jd^x Tneff{x)neff{o) (3) 

Here Tie// is the usual effective hamiltonian describing the low energy weak interactions of 
b and c quarks. The Be state in is to be taken in conventional (relativistic) continuum 
normalization, {Be\Be) = 2EV. 

In the case of heavy quark decay, where the energy release is large, one may perform an 
operator product expansion (OPE) of T. In this way the expression in (^ is expanded in 
a series of local operators of increasing dimension, whose contributions to P are suppressed 
by increasing inverse powers of the heavy quark masses. This formalism has already been 
applied to calculate the total rates for charm and bottom hadrons containing one single 
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heavy quark [0]. For earlier work using similar methods see also [0,0]- Here we shall 
extend this approach to the treatment of B(., where both constituents can be considered as 
heavy. The operator product expansion leads to the generic result 



T = J2CMOn{fi). (4) 

n 

We now describe the operators that contribute to corrections up to order to the free 
quark decay, where v is the quark velocity in the Be meson. For simplicity of discussion we 
will for a moment not distinguish the b and c quark velocities and masses (in the following 
Q can be either b or c). The velocity scaling rules themselves will be discussed in the 
subsequent section. The dominant contribution to the decay width is generated by the 
operators O^q = QQ. To leading order in the velocity expansion of its matrix element, this 
operator reproduces the free quark decay width of the quark Q. Using the equation of motion 
for the heavy quark fields to eliminate redundant operators, no operator of dimension four 
remains. There is one operator of dimension five, Ogq = Qgcr/iuG^^Q, whose contribution is 
suppressed by relative to the leading contribution. The most important correction comes 
from the dimension six operators O4Q = bTccT'b, which scale as (F, F' collectively denote 
the Dirac and color structure). Finally, the dimension six operators Oqiq = Qafj,u'ypD'^G^'^Q 
and Oq2Q = QD^iG^'^T^Q contribute also at order v"^. In the following we do not compute 
the coefficient functions of the latter two operators, so that the expansion will be complete 
to order v^. We include the operators Ogq and take their contribution as indicative of the 
error due to neglect of other contributions of order v^. 

Performing the OPE results in the expression, 

^ = %bb + %5c + %,PI + %,WA, (5) 



%,WA — Tcs + Tud + X! (^) 

I 

where the first two terms in (|) account for the operators O^q and Ogq and the other two 
for the four fermion operators C4Q. Explicitly, 

^356 = Ffe,,p,,66 - % [2Pel + Perl + K^b{Pcl + Peel) + i^2fe(Pc2 + Pee2)\ Og^. (7) 



F, 



%^e = re,.pecCC - ^ [(2 + Koe)Psl + i^2cP.2] ^Gc, (8) 

where 
and 

KoQ = c^ + 2cl K2Q = 2{cl-cl). (10) 
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The phase space factors Pi are given as follows []7|,[l4|: 



(11) 



PcTi = \Jl -2{r + y) + {r - yY 1 - 3(r + y) + 3(r^ + y"^) - - y'^ - Ary + 7ry{r + y) 
1 — r — y + Jl — 2(r + y) + {r — yY 



+ 12r\^ In 



Ary 



(12) 



ccl 



1 - 4?/(l -Qy + 2y^ + 12y^) + 2Ay^ In 



1- v^r^ 



(13) 



P 



cc2 



1 - %(1 + I + 3y2) - 3y(l - 2y') In + '^^^ 



1- v^r^' 



(14) 



where ?/ = ml/ ml and r = m^/ml- P^i (-P52) is identical to P^i {Pc2), except that in this 
case y = m'^/m^. Note also that P^ri = Pci for r = and Pcri = Pcci for ^ = y. In decays 
of the b quark, we neglect terms of order m'^/m'^ and set = 0. Furthermore 



T6,P,= ^|K.|V-(1-^- 

47r 



(C+ + c'^_)ibibi)v-AiCjCj)v-A + (C+ - C^_){bibj)v-A{CjC,i)v-A 



(15) 



■gap + — 



12 



pI 



(c+ - c^Yibibi)y_^{cjCj)1^_^ + (5c+ + c^_ + 6c+c_){bibj)y_^{cjC^)1^_ 



(16) 



%r = -^\V^,\'pI ■ 



n 



(1 - Zr f 

12 



-fl'a/3 



'(l-2;r)^ (1 - Z^)^\ 



pI 



{b 



■i^j)v-A{^j^i)v-A 



(17) 



0) 



p±=Pb± Pc 

In the equations above the QCD correction factors have been written generically 



as(/i) 



(18) 
(19) 



(20) 
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where / is the number of flavors. The scale /i is approximately rrib in and rric in 
respectively. For (15) and (|l^) one might anticipate a scale /i ~ 2mhmc/{mf, + rric), with 
twice the reduced mass representing the characteristic scale of the be bound state. Of course, 
the scales are not flxed precisely and one has to allow for a certain variation in /i, indicating 
an uncertainty that is due to neglected higher order QCD effects. A clariflcation of this 
issue requires the consideration of next-to-leading corrections. 

The contributions of the leading operators bb and cc correspond to the imaginary part 
of the diagrams in Fig. ^ which are contained in expression (^. The coefficients of bb 
and cc in ([7p, can be obtained in the usual way by matching the diagrams of Fig. |I|, 
corresponding to the leading terms of the full expression (|^), onto the operators bb, cc. 
These coefficients are equivalent to the free quark decay rate and are known in next-to- 
leading logarithmic approximation in QCD |l5| -p!9| including the charm quark mass effects 
to 0{as) |jT9|. To include the next-to-leading log effects the Wilson coefficients in the 
effective weak hamiltonians are required at next-to-leading order and single gluon exchange 
corrections to the diagrams in Fig. |l| need to be considered. The complete next-to-leading 
order corrections are incorporated in the numerical analysis below and denoted by Tq^spec 
in d^) and (^. Similarly the contributions with Ogq are obtained when an external gluon 
line is attached in all possible ways to the inner quark lines in Fig. The corresponding 
coefficients are known in leading log approximation. Finally, the dimension six operators 
and their coefficients arise from all those contributions, where one of the internal charm 
lines is 'cut' in the diagrams of Fig. |l]a. The resulting graphs are depicted in Fig. 0. These 
contributions are also known in the literature as 'weak annihilation' (Fig. 0a) and 'Pauli 
interference' (Fig. ^ - note the orientation of the fermion lines). 

The expressions (^, (P) have been derived in 0] (see also [^) and are also discussed 
in [[^. The coefficients for the dimension six operators given in (^)-([T^) are new. They are 
valid to leading logarithmic accuracy in QCD and include charm quark and r lepton mass 
effects. 

In order to obtain the total decay width, the expansion of the transition operator (^ 
has to be taken between the meson states. This step involves the calculation of matrix 
elements of the local operators listed above. In general, this is a difficult non-perturbative 
problem which in practice may introduce considerable uncertainties in the evaluation of the 
hadron lifetime. This is the case for instance for heavy-light 6 or c flavored mesons |]7|,p2|. 
By contrast, an important advantage of the present system consisting of two heavy quarks 
is the applicability of a nonrelativistic treatment. Given the successes of this approach 
in describing cc and bb quarkonia, one can expect a rather reliable determination of the 
required matrix elements using the nonrelativistic formalism, either within the framework 
of potential models or by employing lattice QCD. As we shall see, this treatment transforms 
the l/rriQ expansion into an expansion in powers of the heavy quark velocities. 

Before discussing the evaluation of matrix elements and the numerical results for the 
Be meson lifetime we address briefly the physical interpretation of the various terms in (^ 
- (pi^l). First, let us consider the strictly asymptotic limit in which mq/KQCD ^ oo. In 
this limit the be system becomes extremely nonrelativistic, consisting of two weakly bound 
heavy quarks, slowly moving in a Coulomb type potential. The total decay rate for that 
system is then simply given by the sum Fb^^pec + ^c,spec of the weak decay rates of the 
quasi-freely and independently decaying heavy quarks. In the asymptotic limit we have 



5 



(-Be I I -Be)/ (2Mb J = 1 + 0{vq), Q = b, c, and we indeed recover this simple picture 
within the formahsm described above. Of course, the truly asymptotic case is unrealistic. 
However the simple spectator decay contribution is seen to emerge as the formally leading 
term in the operator product expansion. Bound state corrections are described in by 
the operators Oqq (and others of higher dimension which we omitted) and corrections to 
the asymptotic value of the matrix element of QQ. Since the contributions due to the 
chromomagnetic operators Ogq can be related to the Be — B* mass splitting, it is actually 
of order in the nonrelativistic expansion and thus formally of sub-leading importance. 
The only correction of order v'^ arises from the kinetic energy of the quarks inside the bound 
state, which leads to a reduction of the matrix element of QQ by a factor 1 — Vq/2, obviously 
representing the effect of time dilatation. 

The question of bound-state corrections to the decay rate of a system of two heavy 
fermions has also been studied, from a different point of view (using a Bethe-Salpeter 
formalism), in (see also [Q). It has been shown in this work, that the leading net 
correction to the rate due to bound-state effects is just from time dilatation, — 1>^/2. Other 
conceivable contributions, like phase space suppression and Coulomb enhancement turn out 
to effectively cancel in the final answer. As we have seen, this picture is in accordance with 
the findings obtained here within the OPE approach, which incorporates this result in a 
rather simple and straightforward manner. 



III. MATRIX ELEMENTS 

Let us next complete and summarize the discussion of the required Be matrix elements, 
part of which we have already encountered above. The operators described in the previous 
section are still expressed in terms of four-component fields Q. In a system containing a 
nonrelativistic quark, anti-quarks can not be produced, since this would require an energy 
larger than mq. It is then appropriate to 'integrate out' the small components of the field 
and express the result in terms of a two-spinor ipQ. In this way, all contributions from 
scales larger than /i, where ttlq > fi > mQVq, are made explicit and can be accounted for 
perturbatively. This procedure is standard from nonrelativistic approximations to QED and 
Heavy Quark Effective Theory and leads to 

QQ = i^h^Q - ^^U'Dfi^Q + 

y^h9^ ■ B^Q - T4^i(^ ■ 9E)^Q + ■■■ (21) 



2m^ ^ 4mQ 

Qgcj,.G^''Q = -2ij^Qga ■ Bi^Q - ^i^U^ ■ g^Q + ■■■, (22) 

where we have already omitted the term ipga ■ {gE x D)ipQ (spin-orbit coupling), whose 
matrix element vanishes in the valence Fock state of a pseudoscalar meson. The two-spinor 
ipQ is here defined to have the same normalization as Q, 

I d^x ^l^ijQ = I d^'x Q^Q (23) 
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To the required order i/jq is then related to the upper components of Q 



through 



(24) 



+ (25) 



which can be verified by using the equations of motion. Note also that the covariant deriva- 
tive is understood to be in the adjoint representation when acting on the chromoelectric 
field 

{D-E) = {Bt" - gr^^T^A^)E'' (26) 



Equations (^) and (^) are valid up to terms of 0{v ). In the numerical analysis below we 



shall also neglect the small term of 0{'p^ /itlq), which yields the 0{v^) correction to time 



dilatation 

^ ~ 2m^ ' 8m^ ' ■ ■ ■ ^ 2""^ 8 
Radiative corrections modify the coefficients of the chromomagnetic (ff-B) and the 'Darwin 



{D ■ E) term in (pTI) . In the present context these effects can be neglected consistently. 

The matrix elements of the operators on the right hand side of (^Tj) and (^) can be 
evaluated in the nonrelativistic effective theory for example on the lattice, the obvi- 
ous advantage being the possibility of using coarse lattices, since short-distance effects are 
already accounted for. To obtain a desired accuracy in f^, the appropriate number of cor- 
rections to the leading order effective Lagrangian have to be retained. Preliminary studies 
of B(. mesons on the lattice have recently appeared |^|, but are not yet competitive with 



phenomenological potential models. 

To assess the importance of various contributions, we recall that the quark field ipQ 
scales with the heavy quark three velocity as (mQfg)'^/^, a spatial derivative as uiqVq, the 
electric field gE as friQVQ and the magnetic field gB as toqVq (in Coulomb gauge). The 

1 12 

coupling g m a. matrix element counts as Vq . In the present case of a system with unequal 
quark masses, one has to keep in mind that additional factors from the mass ratio can 
enhance or suppress a given contribution. These scaling rules imply that the last term kept 
in (|21| ) and (|22| ) is of the same order as the chromomagnetic term a ■ B, in contrast to the 
analysis of heavy-light mesons where this term is suppressed by Aqcd/'"^q relative to the 
chromomagnetic interaction. The velocity counting relies on the inequality p ~ 1 GeV > 
^QCD; where p is the typical quark momentum in the bound state. 

Let us now turn to the evaluation of matrix elements in potential models for the bound 
state. Denote by T = mcv'^/2 + mbvl/2 the average total kinetic energy of the quarks in the 
bound state. Then, since mbVb = riicVc, we obtainQ 



*To avoid additional notation, the matrix elements are written for b quarks rather than antiquarks. 
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2Mb, ■ ml mbirric + nib) 



2Mb^ ■ ml nicinic + m,) 



T. (29) 



The kinetic energy T is known rather precisely, within about 10%, T = 0.37 GeV |^. It is 
also approximately independent of the reduced mass of the bound state (for the 'logarithmic 
potential' this statement is exact). This value of T implies v1 ~ 0.38 and vl ~ 0.035. A 
straightforward calculation yields the matrix elements 



2Mb, 3^ nib 



(BMUD ■ qEUJBc) 4 , 

2Mb ^ 1^(0)1 ' (31) 

where ^E'(O) denotes the wavefunction at the origin. The matrix elements for 6-quark fields 
are obtained by interchanging m, and m, on the right hand side. The second equation can 
also be obtained from using the equation of motion for the chromoelectric field. The matrix 
element of the resulting four fermion operator factorizes to leading order in f ^ and can be 
evaluated in a straightforward way. 

We note that in the potential model \I'(0), the decay constant Jb, and the vector- 
pseudoscalar spin-splitting are related by 



fi^iimi M.,-M,^.yi^. (32) 

Mb, 9 iTibnic 

For the numerical analysis to follow we take Mb* — Mb, = 73MeV from With the 
parameters of the Buchmiiller-Tye potential this spin-splitting implies J'b, = 500 MeV, which 
we use as our central value for the decay constant. 

Combining these results, we estimate the matrix elements (pT|) and (0). Notice that 
the matrix element of cga ^^G^^ c is dominated by the divergence of the chromoelectric field 
rather than the spin-spin interaction, because the latter is suppressed by mc/mj,. It is 
interesting to detail the deviations of the matrix element of QQ from unity for the c-quark, 
where relativistic corrections are the largest. We find 

{B,\cc\B,) ^ I . 3 Mb*,- Mb, nib 



2Mb, 2 ' 4 mc \ 2mJ 

^ 1 - 0.190 + 0.037-0.061 + .... (33) 

Despite the large c-quark velocity, this expansion appears quite well-behaved. As antici- 
pated, the largest reduction of the rate is due to time-dilatation, representing the dominant 
bound state effect. Furthermore we have 
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The matrix element of bb and bga ^^G^^b can be obtained by interchanging mi, ^ rric in the 
equations above. 

Finally we need the matrix elements of the four quark operators in (^) - (|17|). To 
estimate their values, we employ factorization. This method, which lacks a firm footing 
for heavy-light mesons, can in fact be justified in the present case of B^. Deviations from 
factorization arise from higher Fock components of the wavefunction and therefore are of 
higher order in the nonrelativistic expansion. The relevant matrix elements are then given 
by 



{B,mb,)'^_Ac,c,)^v^ABc) 



(35) 



{Bc\{bA. 



[CjCj)y_ 



2 a/3 



(36) 



where q is the Be meson four- momentum, q = M^^ ^ (m^ + rric) , and fs^ is the Be meson 
decay constant (in the normalization in which /^r = 131 MeV). Deviations from factorization 
modify the decay rate only at order relative to the free quark decay. Using (|35D and (|36[) 
the matrix elements of the dimension six contributions to the transition operator are found 
to be 



{Bc\%,pi\Bc) 
2Mr„ 



12^ 



2ci 



(37) 



{Bc\Tcs\Bc) 

2Mn 



2471 ' •'^c " '^^ 



4cl + cl+ 4c , 



(38) 



{Bc\%t\Bc) 



2M, 



Be 



<,7T 



(39) 



The matrix elements of Tud, T^e and 7^^ are negligibly small as a consequence of helicity 
suppression. 

In the following section, we evaluate the QCD correction factors c+/_ in the above 
expressions at a scale /i = 2m^ed ~ 2.3 GeV, that is characteristic for a transition involving 
both bound state quarks. If this scale were widely separated from m;,, one would switch 
to the effective theory at a scale /i ^ m^, and scale the operators with their anomalous 
dimension to the lower scale. Since 2mred ~ ^b/2, we excercise our freedom to choose the 
matching scale within some variation of rrib and evaluate the QCD correction factors directly 
at 2?77'y.g(^. 
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IV. DISCUSSION 



We are now ready to collect the various contributions presented above and to derive an 
estimate for the meson lifetime. We first summarize the input parameters that we will 
use in our numerical analysis: 

rUb = 5.0 GeV = 1.5 GeV = 0.2 GeV iKbl = 0.04 (40) 



Mb, = 6.26 GeV M^, - Mb, = 0.073 GeV T = 0.37 GeV Jb, = 0.5 GeV (41) 

The parameters Mb,, Mb* — Mb, and Jb, are taken from and T is taken from 0. The 
numbers for mf, and rric correspond to pole masses, the figures in (|40D representing our 
central values. We comment on their choice below. The renormalization scales fi are chosen 
as follows: fii = nib in decays of the b constituent, /i2 = rric in decays of the c constituent 
and /is = 2mred in modes involving both bound state quarks. 

As already mentioned, the pure spectator decay rates Tb^spec and F^spec are known to next- 
to- leading order in perturbative QCD [1T5|-[I9|1. The most complete calculation, including final 



state mass effects in the QCD corrections, is due to [|T^]. We use their results to evaluate 
the spectator quark decay contributions. Thereby we take proper account of charm quark 
and tau lepton mass effects. The strange quark mass is kept for the charm decay modes, 
but neglected in 6-decays. The impact of electron and muon masses is likewise small and 
has been neglected throughout. 

First we would like to present an overview over the different Be decay mechanisms and 
their relative importance as obtained within the nonrelativistic OPE based framework we 
are advocating. To this end we fix all input parameters at their central values and calculate 
the partial, semi-inclusive Be decay modes corresponding to the various underlying quark 
subprocesses. The results are collected in Table I. We observe a dominance of the charm 
decay modes over 6-quark decay. Weak annihilation is sizable, yet still considerably smaller 
than 6-decay. Adding up all contributions yields a total width of F^^ = 1.914 ps~\ implying 

rB^ = 0.52ps (42) 

Various branching fractions can also be inferred from Table I. For instance the semileptonic 
branching ratio B{Bc — > Xeu) is found to be about 12%. 

These estimates involve considerable uncertainties and Table I is merely intended to con- 
vey the general trend and typical numbers for our favorite parameter set. The dominant 
uncertainty comes from the quark masses. To limit the related ambiguity in our phenomeno- 
logical analysis we employ the following strategy. We take the charm quark mass rric as the 
basic input parameter, allowing it to vary within 

1.4 GeV < me < 1.6 GeV (43) 

which is in the ball park of the pole mass values available in the literature. For any given 
value of rric, "we then fix the value of rrib by the requirement that the measured B^ meson 
lifetime tb^ ~ 1.55 ps is obtained. This is justified since the total B^ width is essentially 
determined by the 6-quark decay contribution with very small pre-asymptotic bound state 
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corrections and the OPE formalism can be expected to be reliable in this case. For the 
6-quark spectator decay rate we again use the complete next-to-leading order expressions 
|T9| . For completeness we also incorporated the bound state effects that are described in [|^ . 
Imposing the constraint yields a 6-quark mass that effectively includes unknown higher 
order perturbative corrections in the decay of the bottom quark. It practically eliminates 
the nib and Vcb dependence of the predicted value of tb^. It turns out that a determination 
of rrib, for any given rric, via the approach just described, is approximately equivalent to the 
relation 



mb = mc + 3.5 GeV, (44) 

which is roughly consistent with the well known formula relating nib and rric in HQET. 
Adopting this procedure and varying nic between 1.4 and 1.6 GeV results in the prediction 

0.4ps < < 0.7ps. (45) 



Variations of our central value ( ^21) due to variations of other input parameters can be 
inferred from Table This table could also be used, if by the time tb^ is measured a 
different parameter set appears preferred. 

In addition to using the measured Bd lifetime one could be tempted to use the lifetime 
to eliminate the strong dependence on m^. Since a very large charm mass is required to 
reproduce the absolute lifetime and semileptonic width in the OPE approach, this would 
lead to the prediction of a very low lifetime tb^ ~ 0.35 ps. We refrain from this procedure, 
since the OPE applied to charmed mesons is less reliable than for Be and probably more 



qualitative than quantitative |22 



Let us conclude the discussion with the following remarks. 

(a) As already discussed above, the expression for the Be meson total width we have 
derived is consistent up to and including terms of 0{v^) in the nonrelativistic expansion. 
Part of the 0{v'^) corrections, those due to the operators Qgcr^^G^'^Q, have also been con- 
sidered to obtain an estimate of the order of magnitude of these contributions, while the 
neglected contributions arise from the operators OeiQ and Oq2Q listed earlier. The labor 
involved in calculating these remaining contributions does not appear justified in view of 
the uncertainties connected with the quark masses. To reinforce this point and to see the 
behaviour of the expansion more explicitly, we write down the size of the various terms 
of a given order in f , using central parameter values. For the sum of the c — > s decay 
modes one obtains (1.526 — 0.289 — 0.008) ps~^, writing separately the terms of order 

and v^. For the total 6-decay contribution one finds (0.640 — 0.011 — 0.014) ps~^. The 
sum of all effects amounts to typically 0.1 ps~^. In the case of 6-decay the corrections 
are particularly small, since they are additionally suppressed by inverse powers of the large 
6-quark mass. Generally speaking, the relevant velocity is the one of the lighter constituent, 

which determines the convergence properties of the nonrelativistic expansion. As we see, 
the series is rather well behaved. Note also that a further contribution of 0{v'^) comes from 
the expansion of the time dilatation factor. Numerically, — ~ —2%. Although one has 
to be careful about drawing definitive conclusions, since the 0{v'^) contribution has not yet 
been calculated completely, these observations support the assumption that the total 0{v^) 
term is rather small. 
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Let us re-emphasize that, technically, in the approach adopted here the velocity expansion 
distinguishes a heavy-heavy meson from a heavy-light meson, for which an expansion in 
Aqcd/'"^q is appropriate. This distinction was missed in the treatment of I^G] (which also 
starts with an OPE), resulting in an incorrect evaluation of the matrix elements of Qga-GQ. 
As a consequence a large correction was obtained from this formally subleading {0{v'^)) 
contribution, which we found to be essentially negligible, even for charm. The leading 
correction terms to the spectator picture of 0{v^) from the kinetic energy and of 0{v^), 
from weak annihilation and Pauli interference, were not calculated explicitly in ||26|| . 

(b) The Pauli interference contribution exhibits a fairly substantial dependence on the 
renormalization scale. The number shown in Table |I| corresponds to fi = 2mred ~ 2.3 GeV. 
Allowing a variation of /i from 1 — 5 GeV yields a range of values from — 0.342 ps^^ to 
— 0.036 ps~^ for this contribution. This large dependence is formally of 0{as), which goes 
beyond the leading log approximation we are working in presently. It represents a theo- 
retical uncertainty in this calculation. The pronounced sensitivity to the scale arises from 
sizable cancellations that occur between the Wilson coefficients in (|37|). Note that the Pauli 
interference contribution is positive in the limit ^ and changes sign due to the pres- 
ence of important short- distance QCD effects. The situation could be improved by studying 
next-to-leading order QCD corrections to Pauli interference, including a proper matching of 
the operators to non-relativistic QCD. By contrast, the scale dependence turns out to be 
very moderate in the case of weak annihilation of be into cs, only ±10% for the same range 
of /i as before. 

(c) For the channels b — >■ ccs and b — > crz/ the weak annihilation term can be comparable 
to, or even exceed the spectator contribution. This feature is related to the strong phase 
space suppression of the three body decay mechanism, which is absent for weak annihila- 
tion. It would persist even in the limiting case where is taken to be very small and the 
nonrelativistic bound state description would be a perfect approximation. The charm mass 
could be so large that the spectator decay would be kinematically forbidden in which case 
weak annihilation would just be the leading contribution, but still reliably calculable in the 
nonrelativistic approach. For these reasons the dominance of weak annihilation does not 
indicate a problem for the validity of the operator product expansion. 

(d) An important conceptual point that needs to be mentioned is that the application of 
the OPE is based on the assumption of local quark hadron duality. By this we mean that 
the sum over all decay channels can be described in terms of partonic degrees of freedom. 
Although this assumption should be valid asymptotically as mg/AQCD oo, it is not a 
priori clear how well it is satisfied in practice. Little is known rigorously and quantitatively 
about the exact conditions to be met if duality is supposed to hold. Conceivably, the issue 
of duality is related to the fact that the velocity expansion is at best asymptotic. This 
includes the possibility that it might not even be asymptotic at minkowskian kinematics. 
Since, for B^, the expansion is still convergent, we are concerned mainly with the second 
possibility, in which case the OPE is meaningless beyond a certain order, even if it appears 
to be convergent. Although duality has been questioned for the decay of D mesons it 
seems, at least to the authors, that there is as yet no place in the charm and bottom hadron 
family where this possibility could be unambiguously separated from uncertainties due to 
unknown input parameters (quark masses, matrix elements) or convergence properties of 
the expansion. In this situation one has to rely on a more pragmatic attitude towards 
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the problem, assume that quark-hadron duahty makes sense for a given case, explore the 
consequences and, if possible, check whether the emerging picture is at least consistent. 

The critical point with regard to this issue for is charm decay. Here the energy release 
is not as comfortably large as it is in the case of bottom decay, and in fact slightly smaller 
than in D meson decay. However, the reliability of a duality is channel dependent. Duality 
might still work reasonably well in the case of the charm transitions contributing to Be 
decay, even if it were violated for D mesons. After all, as we have emphasized, the nature 
of the OPE is very different for Be due to the nonrelativistic character of the heavy-heavy 
type bound state. 

A useful cross-check comes from comparing the inclusive decay width for a given channel 
with the sum of the corresponding exclusive decay modes. In the present case, c sud, 
these are modes like Be B^*^7i, B^*^p, B^*^K^*\ .... In |QJ^ phenomenological models have 
been used to estimate the rates for twenty of these two-body decay modes. Adding their 
results one finds a total of 0.49 ps~^ (BSW model) and 0.65 ps~^ (ISGW model) in |^ and 
0.67 ps"^ in [H, to be compared with 0.91 ps^^ for the inclusive width from Table |. Viewing 
this comparison with due caution, regarding the model dependences and other uncertainties 
in the estimation of exclusive modes as well as the charm mass uncertainty for the inclusive 
prediction, it is reassuring that, first, the order of magnitude comes out to be consistent 
and, second, the sum of exclusive modes does not exceed the inclusive result. A similar 
observation holds for the semileptonic c — > seu transitions. 

Although the arguments we have given are of a somewhat heuristic nature, they all seem 
to indicate that the OPE approach makes sense for inclusive Be decays. At the very least 
the underlying assumptions do not seem to be obviously violated. 

(e) In the Be lifetime has been estimated on the basis of a modified spectator model, 
where the phase space for the free quark decay is modified to account for the physically 
accessible kinematic region or the b and c quark masses are reduced by the binding 
energy to incorporate bound state effects j^. Like the OPE formalism these approaches 
lead to a reduction of the free quark decay rates caused by binding, but the details of the 
resulting picture are markedly different. The modified spectator ansatz does not correctly 
approach the asymptotic limit of very heavy quark masses, where, as in the OPE-based 
calculation, the leading bound state corrections come only from time dilatation. Our bound 
state corrections are numerically much smaller than the very large effects reported in 0, 
which invert the hierarchy between c and b decays, making the latter more prominent, and 
lead to a considerably longer lifetime of tb^ = 1.35±0.15 ps [Q. It has been emphasized in 
that the model calculation has to be regarded with caution, since the inclusive charm decay 
rate is found to be below the rate estimates of the corresponding exclusive modes. Due to 
this feature the result for the lifetime quoted above has been corrected to tb^ = (1.1 — 1.2) ps 
This estimate is however still sizably larger than the one derived in the present analysis, 
where, as we have seen, the inclusive and the exclusive approach can be viewed as roughly 
consistent. 

V. SUMMARY 

In this article we have presented a detailed investigation of the Be meson total width 
based on a systematic operator product expansion of the transition operator. In several im- 
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portant aspects this analysis goes beyond the estimates derived previously in the literature: 
First, we have emphasized the nonrelativistic nature of the heavy-heavy bound state system 
Be- As we have shown, this fact has crucial implications for the organization of the operator 
product expansion, which is different from the case of heavy-light mesons. In addition the 
framework allows a systematic evaluation of the relevant matrix elements of local operators 
using the nonrelativistic bound state wave functions. Matrix elements of four-quark opera- 
tors, for instance, factorize exactly to leading order in the velocity expansion in contrast to 
the general situation. 

Pree-quark decay of h and c represents, both formally and numerically, the leading (in 
v) contribution to the decay rate. For their evaluation the complete next-to-leading 
order expressions in renormalization group improved QCD perturbation theory have been 
employed. The dominant bound state corrections come from time dilatation and are de- 
termined by the heavy quark velocities, which we have estimated from existing potential 
model calculations. A further new result are the complete 0{v'^) contributions to the decay 
rate, due to weak annihilation and Pauli interference, for which explicit expressions, valid 
to leading logarithmic accuracy in QCD, have been derived. 

Numerically our analysis yields tb^ = 0.4 — 0.7 ps, where the measured Bd lifetime has 
been used to reduce the uncertainties due to the bottom quark mass and V^^. 

We have also briefly considered the reliability of the operator product expansion, whose 
applicability relies on the assumption of quark-hadron duality. Although this issue remains 
an important caveat to be kept in mind, we have found no indication for an obvious violation 
of its validity. On the contrary, the series itself is rather well behaved and model calculations 
of exclusive modes are consistent with the duality assumption. It will be interesting to 
compare the predictions for B,, decay based on local duality with future experimental results. 
These should help to assess to what extent the underlying theoretical assumptions are valid 
in practice. Clearly, experimental progress towards a measurement of the lifetime and other 
inclusive properties of Be is eagerly awaited. 
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Mode 


Partial Rate in ps ^ 


b — > cud 


0.310 


b — > CCS 


0.137 


b —f a u 


0.075 


b — > criy 


0.018 




0.615 


c sud 


0.905 


c — > sev 


0.162 


J2c^ s 


1.229 


WA: be cs 


0.138 


WA: be TV 


0.056 


PI 


-0.124 


Total 


1.914 



TABLE I. Contributions to the meson decay rate in ps~ . The partial h and c quark decay 
rates are understood to include the nonperturbative corrections from dimension five operators 
as discussed in the text. Cabibbo suppressed modes are implicitly taken into account with the 
corresponding allowed channels. The input parameters are specified at the beginning of Sect. IV. 
WA=weak annihilation, PI=Pauli interference. 



X 




Ix 


rrib 


5.0 GeV 


-1.96 


rric 


1.5 GeV 


-3.10 


rUs 


0.2 GeV 


0.14 




0.04 


-0.72 




5.0 GcV 


0.02 




1.5 GeV 


0.13 




2.3 GeV 


-0.07 


T 


0.37 GeV 


0.16 


Mb* - Mb, 


0.073 GeV 


0.01 




0.5 GeV 


-0.07 


Mb, 


6.26 GeV 


-0.04 



TABLE II. Parameter dependence of tb, against (small) variations of the parameter set {Xc} 
specified at the beginning of Sect. IV and repeated in the second column of the table. Ix is 
defined by dTB,{X,)/TB,{X,) = lx{5X/X^), where 5tbAX^) = tb,{Xc + ^X) - tb,{Xc)- Thus, 
for instance, the entry for X = rxic tells us that increasing rUc from 1.5 GeV to 1.6 GeV decreases 
tbX^c) = 0.52 ps by 20.7%. 
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Fig. 'Lb 

FIG. 1. Leading contributions in to the spectator decays. 




FIG. 2. 'Weak annihilation' (a) and 'Pauh interference' (b) contributions to the coefficient 
functions of four fermion operators in the OPE. 
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